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ABSTRACT 


mae theoretical acoustic behavior of an underwater bubble 
screen was examined using Rayleigh reflection theory. A 
microcomputer model simulates the acoustic impedance mismatch 
at the interfaces of an ideal, bubble screen. A sensitivity 
analysis indicates that the angle of incidence of sound 
energy and the speed of sound in the layer are the most 
important screen properties for predicting the insulating 
mea ollity of a bubble screen. In the neighborhood of 
frequencies for which the screen thickness is an integral 
number of half wavelengths, the interference results in a 
reduced reflection coefficient and a corresponding increase 
in transmission through the screen. SO that for a broad 
band spectrum wide enough to cover a number of such 
frequencies, the attenuation to be expected exceeds 10 dB 
only over about 90 percent of the spectrum. The interest 
for this work came from the need for a noise insulating 


screen at the Carr Inlet Acoustic Range. 
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Ls ENE Re Due LLON 


Acoustic trials of ships and submarines are conducted by 
Puget Sound Naval Shipyard (PSNS) at the Carr Inlet Acoustic 
Range in the state of Washington. According to Mr. John 
Kriebel of PSNS (personal communication) the measurements 
are often contaminated by noise generated by ship traffic 
in adjacent waters. One method which may offer noise isola- 
tion is a bubble screen located between the range and the 
offending noise source. Mr, Kriebel suggested that a 
feasibility study be conducted. 

This thesis is one part of the study and has the 
objective of examining the theoretical benefits of a bubble 
screen, using a computer model based on Rayleigh reflection 
theory for an homogeneous layer. Recommendations and con- 
clusions are made based on a literature search anda 
sensitivity analysis of the computer model. A companion 
thesis by LT Kelley examines the noise generated by the 
bubbles themselves. 

The writer wishes to acknowledge the interest and 
Support from the PSNS Carr Inlet Acoustic Range during 


moms Study. 
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ie BACKGROUND 


A. ReEoONANCE BUBBLE THEORY 

Experimental studies of the acoustic behavior of bubbles 
in water have been conducted for many years. Normally, when 
a Single bubble is being studied, the resonance theory of gas 
bubble pulsations proposed by Minnaert in 1933 [Ref. 1] is 
considered. This theory shows that a gas bubble in water is 
capable of vibrations which become large at an exciting 


frequency, £, in hertz given by 


pore ile Bie 


ss 2r . 


where, 


f 


. density of water medium in kilograms per cubic 
meter 


r 


ratio of the specific heats of the gas 


ambient hydrostatic pressure in pascals 


radius of the bubbles in meters 
Table I shows the radius of a resonant bubble as a 
menmetion of frequency £f and depth z of concern at Carr Inlet. 


TABLE I 


Resonant Bubble Sizes 







50 m 100 m 50 m 200 


= 








moO Hz Ors 2 mM QO.79 m LOS am oes 

100 O75 03.2. Mm 0.079 m 0.108 m O13) sm O25. m 
LOO0O 3.2 mm 7.9 mm LO o mm 13 mm 15 mm 
10,000 0.32 mm 0.79 mm 1.08 mm 1.3 mm 1.5 mm 


wioeseta 1.93 x hon maya 
I~ 10 Gite O. Lz) pascal 


‘ = 1.40 
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Carstenson and Foldy [Ref. 2] examined certain acoustic 
properties of bubble screens and demonstrated reasonable 
experimental agreement with the theory of resonance of 
bubbles in water. 

Although resonance theory predicts very large attenua- 
tions of incident sound energy for bubbles in water, this 
attenuation is a very strong function of the bubble size, 
"75 percent of the energy lost by the sound wave is due to 
bubbles of radii within 10 percent of the resonant radii." 
[Ref. 3: p.13] A parameter used to describe bubble behavior 
Meets extinction cross section (@%), the ratio of the 
energy lost from the sound wave to the intensity incident on 
the bubble. [Ref. 4] Gis made UD OE DOENHe absorption cross 


2 


section (G ) and scattering cross section (Je) where 


EEG 
For a plane wave propagating ina bubbly mixture, the 


Mmetenvation, a, for bubbles of uniform radius r is 


Mere 5: p.615] 


Q 434-95 -n(r) in JB per meter 


number of bubbles per cubic meter of 
radius r 


where n(r) 


Medwin [Ref. 6] has indicated that "...at a resonance, the 
Betecering and absorption cross sections of a typical bubble 
at sea are of the order 1000 times its geometric cross 


Section.” 
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Several points should be made here. 

meee ne Sensitivity Of attenuation on bubble sizes 
requires a precise control of bubble size. 

2. As will be discussed later, even if bubbles of 
m@aius greater than 1 to 2 cm could be produced, these 
bubbles will rapidly become unstable and break up. 

3. A bubble will expand monotonically as it rises to 
accomodate the decrease in hydrostatic pressure. A bubble 
of 1 mm radius at 100 m will expand to 3.2 mm at the surface 
and its resonance frequency will vary from 10 kHz to 1.0 
mz during its ascent. 

Mec emMenson | Ret) 7) hindicates that 1t is difficult 
to produce bubbles of a given size. Consequently, it is 
mrekicult to predict the attenuation obtainable from given 


mubble production mechanisms. 


i. Bey LEIGH REFLECTION THEORY 

For the above reasons and because of the difficulty of 
producing bubble resonance at low frequencies, this paper 
proposes to view the bubble screen as an homogeneous layer 
separating semi-infinite layers of bubble free water. 
Mmeyreigh reflection theory predicts that attenuation of 
sound energy by the layer of bubbles will be a function of 
Tayer width, angle of incidence, frequency, and the acoustic 
properties of the bubble layer. A layer of constant width 


Will be used to simulate the bubble screen. 
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Peecenslervity sanalysis will then indicate which 
Pemeameters are of greatest concern in bubble layer design. 
Hae sOund power transmission coefficient,%,, will be the 
measure used to analyze these parameters. Lis related to 


the sound power reflection coefficient OX, by the equation, 


ern at Xp =e) 


. MANUFACTURING BUBBLES 

In the past, researchers used a variety of techniques to 
control the size of bubbles in screens. The methods attempted 
Baclude electrolysis, |Ref. 8] microdispersers, [Ref. 9] 
Memos DpOrcelain filters, [Ref. 10] and perforated pipes. 
Exact bubble size was important to determine acoustic 
MmmeaGcteristics of scattering cross section, extinction cross 
section, resonance frequency, and other parameters. Several 
papers concerning bubble generation [Refs. ives; La; 14] 
Mmaatcate that some of the many factors controlling bubble 
mame £Or perforated pipes include: 

ie sOrificie size 

2. pressure drop across the orifice 

meee wetting of the orifice 

4, contaminants in the media 

Bie angeemot Ineclinatilonmof the orifice 

oe rate @£ bubble formation 

Uncertainties in these factors are the primary factors 
fee will prevent accurate bubble size prediction for a 
bubble screen installation on a large scale such as at Carr Inlet. 


a5 





D. BUBBLE HYDRODYNAMICS 

Air bubbles in water exhibit complex hydrodynamic 
characteristics. However, their behavior can be categorized 
generally with respect to the bubble size. Initial bubble 
motion has been analyzed by Walters and Davidson 
[Refs. 15 and 16] and Bachhuber and Sanford [Ref. 17]. 
Gas diffusion from bubble affects bubble behavior signif- 
icantly when the bubble diameter is small. The gas Gilitrusion 
effect depends on the rate of ascent, the radius of the 
bubble and the difference between the gas concentrations 
Meross the bubble Liquid interface. [Ref. 18: ie yaaa 
Mmere are two critical radii: 


Re = the radius above which the bubble grows monotonically 
in size as it ascends; gas leakage is not predominant 


R = the radius below which the bubbles shrink and 
collapse bubbles whose radii lie between Ra and 
Rb are unpredictable. 

Two bubbles of the same gas and size produced at the same 
depth can be affected differently by surfactants, causing 
one of the bubbles to grow and the other to shrink and 
collapse. 

Bubbles of interest for a screen at Carr Inlet are in 
the 1-19 mm diameter range. This size has proven to be the 
easiest to produce via a perforated pipe. This method is 
also by far the most convenient and cost effective. 


"Hydrostatic head has very little influence on the 


memation between the rate of formation and the size of the 
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Supple.” ([{Ref. 19: p.17] VanKrevelen and Hoftijzer 

[Ref. 20: p.30] showed that "...the diameter of the bubbles 
is independent of flow rate and increases with the cube root 
of the orifice diameter." This statement breaks down at 
higher flow rates when gas jetting effects begin. 

Crump [Ref. 21] suggests a differential pressure of 2-3 
psi be established between pipe line pressure and hydrostatic 
water pressure. This should make the desired bubbles for 
forming the bubble screen. A bubble of radius lmm will 
expand monotonically due to decreased hydrostatic pressure 
and will not be affected to a significant amount by gas 
diffusion and temperature effects. A 1mm diameter bubble 
at 150 m depth will increase in size to 2.5 mm in ascending 


l1sothermally to the surface. 


A sample calculation of this effect is shown in Table 
II and shows that this bubble will shift its resonant 
feequency by several kilohertz during its ascent. 


TABLE IT 


Tyowcalbweubble Size Expansion 


RADIUS in mm RADIUS in mm RESONANCE FREQUENCY in KHZ 
mero ft.) Colmer ae eG) (600- ft.) (surface) 

200 a5 TOG Lo. 0 

.500 1.23 agg 6.4 

iF. 0 2.5 ever a2 

Zz. 0 aS tO) iso 

HO . 0 24.6 ise os 
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i. boon OF BUBBLES 
The rate of rise of bubbles has been thoroughly 
Mecumented. (Refs. 22, 23, 24, 25] The speed at which a 
bubble rises depends primarily upon the size and shape of 
mm@enbubble, and can be categorized as in Table Iil. 
Meer. 26: p. 72] 
PAB ee, 2 Lt 


Chasacteristicseor Rising Bubbles 


StAPE DIAMETER in MM SERED. in em s Re NUMBER MOTION 

epherical oo a= © < 200 rectilinear 

ellipsoidal 1-15 10-30 200<Re<4700 Zig-Zag 
helical 


ose Plating 


spherical- mae 30-40 > 4700 errat ike 
capped unstable 


Increasing bubble size causes: 

ale imemeased leecerficient of drag 

Die increased vertical speed 

SeeenonlLinear motion 

ee oscillating motion 

>. wake shedding 
Turbulence seems to be the controlling factor in the motion 
Seebubbles larger than 1200 mm. [Ref. 27: p. 121] "Bubbles 
smaller than about 0.1 cm in water rise along smooth paths 
and do not generate any measureable sound when maser. 


[Ref. 28: p. 25] 
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purtace Contaminants affect bubble behavior in many ways. 
Mhese surfactants [Ref. 29] will have their greatest effect 
on the terminal speed of a bubble by affecting transition 
maemerectilinear motion to oscillating motion. The bubble's 
wake structure is changed causing delayed boundary layer 
separation, delayed vortex shedding and smaller bubble wakes. 

Terminal speed will be a function of bubble diameter, 
surfactants and Reynold's number. Approximate increase in 
vertical speed due to multiple bubbles range up to a factor 


meet. 5 times the single bubble velocities. [Ref. 30] 


_— BuBSpiE NOISE 
Previous investigations [Ref. 31] show that greater than 
PO percent of sound produced by bubbles in water occurs during 


mupble formation. 


Mme © Ff co SC PC. 


10 ms 
migure 1. Sound produced by a bubble at an orifice 
mers sketch [Ref. 32: p. 23] correlates the process of an 
individual bubble leaving an orifice with the observed time 
variation of the radiant noise. The peak sound pressure 
occurs at the time the bubble separates from the orifice. 


Strasberg [Ref. 33: p.24] predicted the peak sound 


pressure at distance d to be, 


i 





_ 1 (gcePvd, |” 
me * er 


where, d = EOtal dnssitpacion constant 
d. = .014 for an air bubble in water 
V = total volume rate of bubble formation 
p = excess pressure 


His works out to be, 
aso dis) ire: 4 ka at tm. 
for a bubble rate of 10/sec for f = 1.0 KHz 

After separation, larger bubbles generate continuous 
sound levels of relatively small amplitude (46 to 66 dB re 
l mPa at lm). [Ref. 34: p.25]| 
This is apparently due to oscillations induced by their 
meregular flow path. Smaller bubbles produce less or 
immeasurable sound when rising. The coalescing of bubbles 


Cause one tenth the sound of bubble formation. 


G. PROPERTIES OF BUBBLE SCREENS 

Studies on the acoustic theory of bubbles and bubble 
Beeeens were done during World War II. [Refs. 35, 36, 37] 
Efforts in the United States were concentrated through the 
National Defense Research Committee (NDRC) under a broad 
research program supported by the Office of Scientific 
Research and Development (OSRD). The prominent NDRC 


Organizations whose findings were published by the U.S. 
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Navy after the war as NDRC Summary Technical Reports were 
mae University of California Division of War Research 
(UCDWR), the Columbia University Division of War Research, 
the Columbia University Sonar Analysis Group, and the Woods 
Hole Oceanographic Institution. The UCDWR was the largest 
of the associated groups and came to be known as the Sonar 
Mme OiVision, accounting for a significant amount of 
empirical studies on ship's wakes, scattering and absorption 
of bubbles, and sound transmission through bubble screens. 
Carstenson and Foldy were part of this effort and in 1947 
they published their study of the propagation of sound 
through a bubble screen. Their screen was approximately 
Meet. Long and 10 ft. high with thickness varying from 

3 in. to 6 in. The pulsed screen was analyzed as if 
bubbles were uniformly distributed between two infinite 
parallel planes and the speed of sound in the screen was 
assumed to be approximately the same as that of water. 

For the case of continuous - flow screens, the bubbles 
were of varying sizes and speed of sound in the layer 
considerably different from that in water. 

Their results supported the resonant bubble theory for 
pulse screens at frequencies from 15 kHz to 35 kHz. However, 
Seem data from 5 kHz to 45 kHz for continuous flow screens 


Of bubbles of varying sizes were not very satisfactory. 


Za 





Throughout their study, the thickness of the screen was 
assumed to be the arithmetic average of the observed layer 
width. The density of the screen was based on observations 
of the average number of bubbles present per unit volume 
and of the average bubble size. 


For the continuous flow case, the screen had a "..fairly 
mmerply Gefined core of larger bubbles...but in front and 
back there is a gradually tapering distribution of smaller 
Pibbles.” [Ref. 38: p.501] As a result, the speed of sound 
massa Gradually changing function of the distance normal to 
the layer interface plane causing reduced reflections of 
Sound . 

Meod and Spitzer [Refs. 39,40] proposed theoretical sound 
emeeds which would occur for given concentrations of air by 
volume in water. Figure 2 illustrates those results and 
experimental verification by others. [Refs. 41, 42] For 
concentrations between Mor and ion readily obtainable in 


the laboratory or in the field, sound speeds as low as 


100 = 300 m/sec are found. 


fs APPLICATIONS OF LAYERS 

There are extensive studies of the reflection and trans- 
Mission of energy in layered media. [Refs. 45, 46] Numerous 
applications of the use of layers range from reflection 
reduction for optical lens coatings to absorbent materials 


mememarchitectural acoustics. The ratio of layer thickness 


to the wavelength of the energy in the layer is a very sig- 


nificant parameter due to the standing wave pattern produced. 
ae 
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Figure 2. 
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Sound speed in air-water Mixtures 
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a. Bwooseeeacane at CARR INLET 


ae OCEANOGRAPHIC CHARACTERISTICS 

The chart of Carr Inlet Acoustic Range in Appendix A 
indicates a large, irregularly shaped area surrounding the 
test site. Tite m—sleseo rts snougniy 5 Im long and 3 nm wide. 

It is proposed to consider in this study a bubble screen 
crossing the mouth of the inlet, a distance of approximately 
2900 m. This position places the screen directly between 
commercial boat traffic to the southeast and the test site. 

The water depth at this position varies from approximately 
450 ft. at the relatively level bottom at the middle of the 
inlet, to relatively steep shore gradients approaching Fox 
moband and McNeil Island as shown in Figure 3. 

The frequencies of interest for attenuation will be 
meeadmband (100 Hz to 10 KHz+t). 

Oceanographic studies [Ref. 47] at Carr Inlet indicate 
that the tides and currents which are likely to be encountered 
will be variable but less than one-half kt. There is no 
Significant evidence for underground springs or bottom 
currents which would cause any significant salinity or 


temperature gradients. Seasonal variations in water 


memeperature and salinity occur due to the annual air tem- 


Bemeature cycles and rain runoff. However, wo —Sitgnifticant 
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CARR INLET — SIDE VIEW 


aa 
1,000 
4 screen 
{in =< 1mm 
depth 
fa f.%. 


CARR INLET - FRONTAL VIEW 





{in =1/4 nm 
depth 
in ft 
Figure 3. Position of bubble screen in sensitivity analysis. 
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maupact 1S anticipated on the sound velocity profile. So, 
for the distances encountered in this application, the 
sound speed in water will be assumed to be constant 


mee— 1500 m/s). 


Pepe SOCREEN POSITION 

Considering these oceanographic factors, the bubble 
screen will be assumed to be produced at the bottom of 
mmeeintet and to rise vertically with constant width. An 
estimation for actual layer width spreading will be 


presented in Chapter V. 
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se Sens Gl iyity ANALYSIS 


A. maEORY 

At frequencies removed from the resonant frequency of 
the bubbles, a bubble layer is expected to reflect sound as 
a result of two interrelated effects on the acoustic 
Meeperties of the medium. [Ref. 48: P.925}] First, the 
compressibility of the water is changed with the introduction 
of the bubbles which in turn changes the sound speed in the 
layer. Second, density will be changed. For the bubble con- 
centrations considered in this thesis, this latter change is 
emall and the principal acoustic characteristic of the bubble 
layer will be its sound speed. The sound speed in bubble 
layers was addressed by Wood [Ref. 49] whose work was con- 
marmed experimentally by Silberman. [Ref. 50: p.90] 

In the present work, the acoustic effect of a bubble 
screen will be analyzed in terms of the sound power trans- 
mission coefficient using a computer simulation of the bubble 
screen. As indicated earlier, this computer simulation will 
be based on Rayleigh reflection theory for an homogeneous 
layer of bubbly water with constant thickness. 

meryvleigh reflection theory is based on a solution of a 
Sescem Of acoustic equations, describing the conditions that 


must be satisfied at each boundary of the layer. 
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mmemooundary conditions at the interfaces are continuity 
Setne acoustic pressure and of the normal component of the 
fercicle velocity. 

For oblique incidence and when the media in front and 
behind the layer are the same, the sound power transmission 


Beeecicient CX.reduces to [Ref. 51], 





f 
ace a a2 a_: a (1) 
(cos (bh) aE (m, +mg) sin (b, h)} 
This coefficient represents the ratio of the intensity of 
the transmitted sound to the intensity of the incident sound. 


The geometry of the problem is shown in Figure 4, and 


a measured from normal 


m — ©2S2 COS o n 
Cc, cose pomiieemface, n = 1,2,3 
/ / Po} 


wer. neGraent. layer 


Gig e OS , 2....screen layer 


Soe. eOXLt Layer 


b, = 2mrf cos ©: 
P) C, 


h = layer thickness 


= Bt ee 


mma oy Snell's law: 


Cy SiN oy 
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Figure 4. Bubble layer reflection geometry 
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is Mice note okewictme and acoustic impedance of the 
layer and depends on the angle of incidence and frequency of 
the incident sound. 

The subsequent sensitivity analysis will indicate the 


theoretical performance expected when parameters vary. 


Bi. Samo lLTiVITY ANALYSIS RESULTS 
l. Parameter description 
Data were taken using the options in Table lV, aA 
change in symbols was required for use in the computer model. 
The new symbols used to represent the parameters are also 


Mmested for the options in Table lV. 


TABLE IV 


Analysis Options 


OPTION ORDINATE ABSEISSA CONSTANT PARAMETER 
bl TS 2 Al A 
2 T3 Al L2 A 
3 TS Al A ee. 
- TS [2 A Al 
5 T3 A L2 Al 
6 -3 A al ee 
where, 
73 = sound power transmission coefficient 
L2 = layer width divided by the wavelength in layer 
Al = incident angle 
A = sound speed ratio C 


layer / eer 
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Each of these options will be addressed in the following 
analysis. Values of A and their corresponding fraction 
volume of air bubbles are shown in Table V, 
TABLE V 
Sound speed ratios for given air concentrations 


Proce -omuby Volume Of 


A alr in water 
0.5 2x 107° 
Op4 ees 
3 5 2 1” 
O. 2 2x 10> 
On Lt eon 


These values of A may reasonably be expected ina 
practical bubble screen. 

The value of L2 was obtained by incorporating a normal- 
Mmepeton factor which forced the data to repeat itself after 
L2 reached a value of 1.0. This was done for ease in 
Meet cing and clarity in data presentation. The normalization 


meteor fF! is found by noting that in Eq. (1), 
x= 1, for bhenT nz 0,2, 


for what values of layer width/wavelength (~———] ¢ 


Min layer) . 


Bei. =nt 
Co 


Srl 





het n = l, 


aL h cos 0, = 1 
2 7 ers 
where Ll - 2 and, by Snell's Law, COS ©, ~ yf Asi & 


a 


Therefore, Ll -1/(2 y1- A® sin OQ, and L2 = WEE 


eee The Computer Program 

The WANG SYSTEM 2200S desk top computer was used to 
run the analysis and provide data output. The program, 
"OPT6", was written in BASIC language and provides a defini- 
tion of the variables and symbols at the beginning of the 
Program. 

The plotting option is selected when the program is 
executed. Each option has its own subroutine which 
interrogates the operator for the bubble screen parameters 
required to produce the selected data output. Step sizes 
For all the parameters are pre-selected. The sound power 
transmission coefficient is plotted vertically against the 
changing option parameters. 

S . Oprren. | 

When the sound power transmission coefficient is 
plotted against the normalized layer ratio, T3 cycles be- 
tween maxima and minima. Minimum sound power is transmitted 
through the layer at odd integermultiples of one-quarter 
wavelengths. Figures 5 to 10 show the dependence of T3 on 


MmemOoL a fixed Al for a variety of values of A. 
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A = SOUND SPEED RATIO 


Pee = INCIDENT ANGLE 
L2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION COEFFICIENT 
Figure 5. Sound power transmission coefficient vs. 
layer ratio, parameter A for constant 
Al = 0” 
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Al 
a2 


is 


Eeeagure 6. 


SOUND SPEED RATIO 
INCIDENT ANGLE 


LAYER WIDTH/ WAVELENGTH RATIO 
SOUND POWER TRANSMISSION COEFFICIENT 


Sound power transmission coefficient vs. 
layer ratio, parameter A for constant 
Al = 20 
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T3 10 


BJ 
A=.1 
3 1G 
eZ 
Pee = SOUND SPEED RATIO 
Al = INCIDENT ANGLE 
L2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION Cee riGien | 
meaqure 7. Sound power, transmission coefficient vs. 
layer ratio, parameter A for constant 
Al = 40 








Pee = SOUND SPEED RATIO 
Po = INCIDENT ANGLE 
L2 = LAYER WIDTH/ WAVELENGTH RATIO 


T3 = SOUND POWER TRANSMISSION COEFFICIENT 


Figure 8. Sound power transmission coefficient vs. 
layer ratio, parameter A for constant 
Al = 60 
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megure 9. 


SOUND SPEED RATIO 
INCIDENT ANGLE 
LAYER WIDTH/ WAVELENGTH RATIO 


SOUND POWER TRANSMISSION COEFFICIENT 


Sound power, transmission coefficient vs. 
layer ratio, parameter A for constant 
Al = 80° 








A = SOUND SPEED RATIO 


A1 = INCIDENT ANGLE | 
Mee = LAYER WIDTH/WAVELENGTH RATIO 
mee = SOUND POWER TRANSMISSION COEFFICIENT 
Figure 10. Sound power transmission coefficient vs. 
layer ratio, parameter A for constant 
Al = 85 
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The following observations can be made: 

a. Regairimbess eon angle Of incidence, if the sound 
meeca ratio is less than 0.1, attenuation of at least l10dB 
can be expected over at least the middle 75 percent of the 
layer ratio range. 

Des Regardless of the sound ratio or the angle of 
incidence, there will always be at least 5-19 percent of the 
layer ratio which will have only negligible attenuation. 

meee Option 2 

When the sound power transmission coefficient is 
plotted against the angle of incidence, the maximum trans- 
mission always occurs at normal incidence and the minimum 
transmission at angles approaching On. Figuress)1 to Ic 
are plots with the layer ratio held constant and the sound 
speed ratio varied as a parameter. 

The following observations can be made: 

a. When the sound speed ratio reaches 0.1, at 
least 10dB attenuation can be expected over all angles of 
incidence and over at least 75 percent of the layer ratio 
range. 

b. The only values of sound speed which will give 
a uniform value of attenuation over all angles of incidence 
are speed ratios of less than 0.1. This situation would 
require the bubble screen to maintain a volume fraction 


~2 
of air to water of greater than approximately 10 2 
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309 60° | 90? 
A = SOUND SPEED RATIO Al 
Al = INCIDENT ANGLE 
L2 = LAYER WIDTH/ WAVELENGTH RATIO 


T3 = SOUND POWER TRANSMISSION COEFFICIENT 


Figure ll. Sound power transmission coefficient vs. 
incident angle, parameter A for constant 
L2 = 0 
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BO? 60° 7o° 


Al 
A = SOUND SPEED RATIO 
Be = INCIDENT ANGLE 
ime = 


LAYER WIDTH/WAVELENGTH RATIO 


es SOUND POWER TRANSMISSION COEFFICIENT 


Figure 12. Sound power transmission coefficient vs. 


incident angle, parameter A for constant 
2S ed 
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m= SOUND SPEED RATIO 


Pie = INCIDENT ANGLE 
L2 = LAYER WIDTH/ WAVELENGTH RATIO 
i = 


SOUND POWER TRANSMISSION COEFFICIENT 


Fagure 13. Sound power transmission coefficient vs. 


incident angle, parameter A for constant 
2a 2 
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fee = SOUND SPEED RATIO 


Pee = INCIDENT ANGLE 
L2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION COEFFICIENT 
Mmegure 14. Sound power transmission coefficient vs. 
incident angle, parameter A for constant 
b2= .3 
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Soe 60° 9G: 


Al 
Pe = SOUND SPEED RATIO 
me = INCIDENT ANGLE 
L2 = LAYER WIDTH/WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION COEFFICIENT 
Mmegure 15. Sound power transmission coefficient vs. 
incident angle, parameter A for constant 
L2 = .4 


44 








Be = SOUND SPEED RATIO 


Al = INCIDENT ANGLE 
L2 = LAYER WIDTH/WAVELENGTH RATIO 
m= SOUND POWER TRANSMISSION COEFFICIENT 
Figure 16. Sound power transmission GoetrErciemt, vs. 
incident angle, parameter A for constant 
b2 = 5 
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ee Option 3 

As in the plots Mow OopelLoen 2, the sound power 
transmission coefficient is plotted against the angle of 
incidence. This time, though, the sound speed ratio is 
plotted as a parameter. See Figures 17 to 22. 

The maximum sound transmitted occurs at normal 
incidence and the minimum always occurs as the angle of 
incidence approaches 90°. 

The following observations can be made: 

a. The significance of a change in the sound 
Speed ratio is very clear in this option. rf Figure Is 
(sound speed ratio = 0.1) is compared to Figure 22 (sound 
speed ratio = 0.5), a dramatic decrease of attenuation is 
Seem tO occur. 

b. At low angles of incidence, values expected for 
attenuation are relatively constant. At high angles of 
incidence, attenuation values are changing rapidly, 
indicating less confidence for a given predicted layer width. 

oS Option 4 

This option plots the sound power transmission 
coefficient against the normalized layer ratio. The angle 
of incidence is the parameter which is varied for constant 
mmies of sound speed ratios (Figures 23 to 29). Again it 
7S seen that values for the sound power transmission 


coefficient cycle between maxima and minima. The minimum 


46 





eeee= SOUND SPEED RATIO 
Al = INCIDENT ANGLE 


L2 = LAYER WIDTH/WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION COEFFICIENT 
Figure 17. Sound power transmission coefficient vs. 

incident angle, parameter L2 for constant 

A = .05 
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A = SOUND SPEED RATIO 


Be = INCIDENT ANGLE 
L2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION COEFFICIENT 
mmere 18. Sound power transmission coefficient vs. 
incident angle, parameter L2 for constant 
A= .l 
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A = SOUND SPEED RATIO 


Al = INCIDENT ANGLE 
L2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION COEFFICIENT 
Figure 19. Sound power transmission coefficient vs. 
incident angle, parameter L2 for constant 
A= .2 
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Figure 


30° 60° 90° 


SOUND SPEED RATIO 
INCIDENT ANGLE 


LAYER WIDTH/WAVELENGTH RATIO 
SOUND POWER TRANSMISSION COEFFICIENT 


20. Sound power transmission coefficient vs. 
incident angle, parameter L2 for constant 
A = .3 
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30° 60°? 90° 


A = SOUND SPEED RATIO 
Al INCIDENT ANGLE 


L2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION COEFFICIENT 
Mmeoure 21. Sound power transmission coefficient vs. 
incident angle, parameter L2 for constant 
A= .4 
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Pagure 


2 0* SO. 90° 


SOUND SPEED RATIO 
INCIDENT ANGLE 


LAYER WIDTH/WAVELENGTH RATIO 
SOUND POWER TRANSMISSION COEFFICIENT 


22. Sound power transmission coefficient vs. 
incident angle, parameter L2 for constant 
A= .5 
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A = SOUND SPEED RATIO 
Al = INCIDENT ANGLE 
ue LAYER WIDTH/ WAVELENGTH RATIO 


es SOUND POWER TRANSMISSION COEFFICIENT 


Figure 23. Sound power transmission coefficient vs. 
layer ratio, parameter Al for constant 
Ay 2 2015 
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fee = SOUND SPEED RATIO 
Al = INCIDENT ANGLE 


L2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION COEFFICIENT 
Figure 24. Sound power transmission coefficient vs. 

layer ratio, parameter Al for constant 

os 
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Pee = SOUND SPEED RATIO 
Al = INCIDENT ANGLE 


L2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION COEFFICIENT 
Meere 25. Sound power transmission coefficient vs. 

layer ratio, parameter Al for constant 

A= .2 
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A = SOUND SPEED RATIO 


Pie = INCIDENT ANGLE 
f2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = 


SOUND POWER TRANSMISSION COEFFICIENT 


Figure 26. Sound power transmission SOerr cient 4ivs. 
layer ratio, parameter Al for constant 
A= .3 
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fe 1.0 


L2 

A = SOUND SPEED RATIO 
Al = INCIDENT ANGLE 
L2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION COEFFICIENT 
mmaure 2/7. Sound power transmission coefficient vs. 

layer ratio, parameter Al for constant 

A = .4 
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Figure 


SOUND SPEED RATIO 
INCIDENT ANGLE 
LAYER WIDTH/WAVELENGTH RATIO 


SOUND POWER TRANSMISSION COEFFICIENT 


Zo oOUNnG power transmission coefficient vs. 
layer ratio, parameter Al for constant 
A = .5 
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Peee= SOUND SPEED RATIO 


Al = INCIDENT ANGLE 
L2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION GCEEFIGIENT 
Figure 29. Sound power transmission coefficient vs. 
layer ratio, parameter Al for constant 
A= .9 
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mmdes OCCUr at layer ratios which correspond to odd 
multiples of one-quarter wavelength. 

The following observations can be made: 

a. Changes in attenuation due to variations in 
the angle of incidence are very slight for angles less 
than approximately 45°, but the attenuation in uniformly 
predictable and by far the greatest at angles greater 
than 80°-85°, 

De. Variations in the sound speed ratio cause far 
greater effects on the attenuation when the sound speed 
ratio is greater than 0.2. 

Meee ODtion 5 

The sound power transmission coefficient is plotted 
against the sound speed ratio in this option. The bubble 
layer ratio is kept constant while the angle of incidence 
is varied. A much different view of the data is seen in 
Figures 30 - 34. 

The following observations can be made: 

an The best attenuation occurs at very low sound 
speed ratios as was also seen in options 1 and 2. 

b. There is a significant amount of attenuation 
gained for angles of incidence past 75 

wee Option 6 
In this option, the sound power transmission 


Beerficient is plotted against the sound speed ratio with 


60 


the parameter of layer ratio varied for constant values 
@peangles of incidence in Figures 35 to 38. 

The following observation can be made: 

a. Layer width ratio curves vary only slightly 
indicating that layer width is not an important parameter. 
mmeaddition, there is still approximately 10 percent of 
meme layer ratio range (L2< .05) virtually unaffected by 


sound speed ratio or angle of incidence. 
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A 
A = SOUND SPEED RATIO 
mee = INCIDENT ANGLE 


Ee 
ie 3 


LAYER WIDTH/WAVELENGTH RATIO 
SOUND POWER TRANSMISSION COEFFICIENT 


Figure 30. Sound power transmission coefficient vs. 


sound speed ratio, parameter Al for constant 
2 =. 05 
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Figure 


SOUND SPEED RATIO 
INCIDENT ANGLE 


LAYER WIDTH/WAVELENGTH RATIO 
SOU | OVER TRANSMISSION COEFFICIENT 


PcounGce powers rransmicsion Coefficient vs. 
sound speed ratio, parameter Al for constant 
ie2e= al 
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Pee = SOUND SPEED RATIO 
AY = INCIDENT ANGLE 


L2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION COEFFICIENT 


migure 32. Sound power transmission coefficient vs. 
sound speed ratio, parameter Al for constant 
B2 = <2 
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ba 10 


A 
A = SOUND SPEED RATIO 
Pay = INCIDENT ANGLE 
me = 


LAYER WIDTH/ WAVELENGTH RATIO 


T3 = SOUND POWER TRANSMISSION COEFFICIENT 
Figure 33. Sound power transmission coefficient vs. 
sound speed ratio, parameter Al for constant 
Lb2 = 3 
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A = SOUND SPEED RATIO 
Al = INCIDENT ANGLE 


eZ LAYER WIDTH/ WAVELENGTH RATIO 
nS SOUND POWER TRANSMISSION COEFFICIENT 


i 


Figure 34. Sound power transmission coefficient vs. 


sound speed ratio, parameter Al for constant 
L2 = .4 
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Pee = SOUND SPEED RATIO 
Al = INCIDENT ANGLE 


L2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION COEFFICIENT 
Figure 35. Sound power transmission coefficient vs. 

sound speed ratio, parameter L2 for constant 

Al = 0 
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A = SOUND SPEED RATIO 
Pa = INCIDENT ANGLE 


2 LAYER WIDTH/WAVELENGTH RATIO 
rs SOUND POWER TRANSMISSION COEFFICIENT 


Figure 36. Sound power transmission coefficient vs. 


sound speed ratio, parameter L2 for constant 
Al = 30 
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Figure 


2 1.0 


SOUND SPEED RATIO 
INCIDENT ANGLE 
LAYER WIDTH/ WAVELENGTH RATIO 


SOUND POWER TRANSMISSION COEFFICIENT 


37. Sound power transmission coefficient vs. 
sound speed ratio, parameter L2 for constant 
Ae Ce 
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A = SOUND SPEED RATIO 


Al = INCIDENT ANGLE 
L2 = LAYER WIDTH/ WAVELENGTH RATIO 
T3 = SOUND POWER TRANSMISSION COEFFICIENT 
Mmeagure 38. Sound power transmission coefficient vs. 
sound speed Tatro, parameter L2 for constant 
Al = 85 


co 


Vv. BUBBLE SCREEN APPROXIMATIONS AND RECOMMENDATIONS 


A. EPReEQUENCY BANDS WITH LITTLE ATTENUATION 

Options 1 and 6 shows that there will be at least 10 
percent of the layer ratio range virtually unaffected by 
a bubble screen. So for any given spectrum, there will be 
frequency bands of small attenuation which occur when the 
screen thickness is an integer multiple of one-half 


wavelengths in the layer. 


Bes LAYER SPREADING 

Answers to questions on layer-width spreading were not 
found in the literature. Forcing parameters which are 
present in Carr Inlet which would cause a bubble screen to 
lose its simulated parallel interfaces with the surrounding 
media are bubble interactions, turbulent bubble motion, 
mageernlet currents. In order to understand these effects, 
it might be better to conduct an experimental study at 
Saree inlet and not a theoretical study because predictability 
rapidly disappears due to turbulent effects when bubbles 
feeow larger than 1 to 2 cm. 

However, an estimate of the layer spreading may be 
taken from a linear extrapolation of the observed laboratory 
results of Carstenson and Foldy. [Ref. 52] Their screen 


peteaad up to 6 in. for a rise of 10 ft. Figure 39 shows 
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Figure 39. Bubble screen spreading effect 





approximate layer thicknesses for a linear extrapolation. 

The corresponding screen volume increase in this case is a 

maetraor of 60 from 600 £t. to the surface. At the same time, 

the bubbles formed at 600 ft. are expanding their volume 

due to decreased hydrostatic pressure by a factor of about 

18 atmospheres. SoOetOseenits lNyeoehetical@ease, the fraction 

volume of air to water changes by a factor of about one-third. 

This could cause a change of speed of sound of 200 - 300 m/s. 
When considering an approximation like this, the number 

of unknowns involved indicate the uncertainty inherent in 

projecting any laboratory layer spreading data to the 


Gmmensions of Carr Inlet. 


ee, EXPERIMENTAL SCALING 

Mme tractors to be considered for proper scaling of 
laboratory experiments are: 

ies Pressure change during bubble ascent 

2. Height, length, thickness of the bubble layer 

3. Angle of incidence and frequency of incident sound 

4, Path lengths of sound 

Of these factors, the pressure change is the most 
difficult to scale because volume expansion of the bubble 
would cause a significant change in the speed of sound. 

Om Site experiments at Carr Inlet appear to be the only 
viable method to test the results of this sensitivity 


emalysis. 
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1). RECOMMENDATIONS AND ALTERNATIVES 

The purpose of this thesis was to investigate the 
benefits of a bubble screen which may offer some isolation 
from underwater noise interference at Carr Inlet. Asa 
result of this sensitivity analysis and further investiga- 
tion, the following recommendations are made: 


1. An alternate screen position as shown in Figure 40. 


® TEST SITE 


\\\ ALTERNATE SCREEN 





Figure 40. Alternate bubble screen position 
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a. The advantages of this position are: 

oe ema ae iEraction Of the piping (one estimate 
Mere piping to cross the inlet for a uniform screen required 
either a large number of in line pressure regulators or 
memest LOO km of piping). 

(2) Shallower depth allowing less dispersion of 
bubble layer. 

(2a eo lOmbteacelom Of the High attenuation pre- 
dicted in this sensitivity at steep angles of incidence. 

iD. The disadvantages include: 

(1) Increased noise interference when screen 
mmoraced closer to the test site. 

(2) Low frequency sound still not affected. 

(35s Dittracteton Of Sound around edges reduces 
any attenuation. Further study for any perforated pipe 
bubble screen mechanism should include pipe suspension or 
anchoring devices, pipe cleaning methods and compressor 
requirements. 

2. Conduct an on-site layer spreading experiment. 
a As an alternative to the bubble screen solution, 
investigate the benefits of adaptive beam forming to null 


out bothersome interference directions. 
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nee. CONncEUS TIONS 


Pe. Mieminsulating capability of a bubble screen at Carr 

Inlet is a function of the angle and frequency of incident 
sound energy and the width and speed of sound of the bubble 
layer. Of these parameters, the angle of incidence and the 


speed of sound in the layer are the most important to consider. 


B. There will be approximately 10 percent of the frequency 
spectrum over which virtually no attenuation of sound will 
occur no matter what speed of sound could reasonably be 


acheived in the bubble layer. 


C. The maximum consistent attenuation which could be 
expected is approximately 10 dB at bubble concentrations 
-3 
greater than about 5 x 10 ‘ 
These conclusions are based on the assumption that the 
air-water mixture in the layer is an ideal, homogeneous fluid. 
As a next step in the improvement of this model, the case 


of a lossy medium should be incorporated to include the 


effects of absorption of a bubble screen. 
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APPENDIX A: BUBBEE SCREEN POSITION AT CARR INLET. 
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com; UTER PROGRAM 


Meee 5; PROGRAM NAME......... RO as MA oe) or elisiu dis pus ea) 08 644 s+ 6 ws 
Peete PUR oe... LOTS ThE POWER TRANSMISSION COEFFICIENT 
meee FOR A PLANE WAVE INCIDENT ON A BUBBLE SCREEN LAYER IN 
wero DESTRED OUTPUT DISPLAY IS SELECTED FROM ONE OF 
ies OPTIONS BELOW 


RSS| DyeJe JOIN IED IOOUN ICIS” Se iC Ile su (Se 
emesis "A LTS ACOUSTIC IMPEDANCE RATIO = R2/R1 

meeneMm 3 "AL" iS THE INCIDENT ANGLE MEASURED FROM NORMAL 
Meee. "AZ" IS THE ANGLE OF TRANSMISSION IN THE SCREEN 
fumes CC" U6TSd6 THE SOUND SPEED RATIO=C2/1500="A" APPROX 
Meets C2" 15> THE SOUND SPEED IN THE BUBBLE LAYER 
Meee; "LL" IS THE RATIO, BUBBLE SCREEN WIDTH/WAVELENGTH 


mere 5; 'LI" IS THE NORMALIZING FACTOR FOR L 


meen eM 3; "L2" = L/LIL 
feces | 1S°”)62S:6 THE SOUND POWER TRANSMISSION COEFF (OFFICER) 
Meee s | 'S" 61S THE SELECTED PLOTTING OPTION 


MOREM - ........000-. EMOAMMMNNGMOPTTONS .......-0+0sceess00s 
150 REM ; OPTION 1 IS T3 VS L2 FOR CONST Al AND PARAMETER A 
160 REM ; OPTION 2 IS T3 VS Al FOR CONST L AND PARAMETER A 
170 REM ; OPTION 3 IS T3 VS Al FOR CONST A AND PARAMETER L?2 
- OPTION 4 IS T3 VS L2 FOR CONST A AND PARAMETER Al 
- OPTION § IS T3 VS A FOR CONST L AND PARAMETER Al 
, OPTION 6 IS T3 VS A FOR CONST Al AND PARAMETER L2 
210 INPUT "OPTION NUMBER =",S 


meeetr S = 1 THEN 300 

meeetr S = 2 THEN 570 

Beer S = 3 THEN 840 

250 IF S$ = 4 THEN 1100 

meee S = 5 THEN 1380 

meee; > = 6 THEN 1700 

0 ere eeu C eae eee) tee OUT ee igisi ss oie ae wee wo el 


290 REM ; THIS SUBROUTINE PLOTS T3 VS L2 (Al CONST), PAR A 
300 INPUT "INCIDENT ANGLE IN RADIANS(EA 10 DEG=.174 RAD)=",Al 
310 SELECT PRINT 213(64) 

320 PRINT " INCIDENT ANGLE=",Al | 
330 PRINT " Cem Oy oul? @0=1.0), PARAMETER AC.1-.9)" 
340 SELECT PRINT 413 

350 PLOT (10) (0,-50,"+") 

SeOe PLOT (10) (50, 0, "+") 

MeePEOT (-500, 0,) 

MeOeEOR A = .1 TO .9 STEP .4 

390 H=0: V=0: H1=0: Vi=0 


HO0 A2 = ARCSINCA*SIN(AL)) 

Seem = A*COS(A1)/COS (A2) 

feeeeri2 = (1/A)*COS(A2)/COS(AL) 

meee = 1/(2*SQRC1-A*A*SINCAL)*SINC(A1))) 


meee 12 = 0 TO 1 STEP .02 
Sere? = 2°#PT*L2*L1*COS (A2 ) 
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Piee= COSCEZ)*COSCBZ2 ) 

Bee = SEN(CB2)*SINCB2) 

T3 = 1/C(CDI1+¢(M14M2)*(M14M2)/4"D2)) 
Ps ieee S00). Y = INTCT3*500) 


PLOT((H-H1), (V-V1) ,HEX(FB)) 

Hl =H: Vl=V 

NEXT L2 

emorT (-H,-V,) 

NEXT A 

Beto 210 

PRE 4b. ck eee eee Wink OlUMMNIR OP MEON 2..2-..0.......+-. 
REM ; THIS SUBROUTINE PLOTS T3 VS Al (L CONST), PARAM. A 
INPUT "BUBBLE LAYER/WAVELENGTH=",L 

SECT PRINT 213(64) 

PRINT " BUBBLE LAYER/WAVELENGTH=",L 

PRINT " TuemCOeNOmvcm rr CO=e9 DEG), PARAM AC.1-.9) 
PEEECT PRINT 413 

por (10) (0,-50,"+") 


PLOT (10) (50,0,"+") 

Ewom ¢-500,0,) 

FOR A = .1 to .9 STEP .4 
H=0: Vs0: H1l=0: V1=0 

FOR Al = 0 TO 1.56 STEP .04 


Pee] ARCSINCA*SINCAL)) 

Deee= A* COSCA1L)/COS(A2) 

fee= (1/A)*COS(CA2)/COSCAL) 

Be = 2*#PLI*L*COS(A2) 

Mee= COS(B2)*COS(B2) 

mee= SLNCB2)*SINCB2) 

jeoe= 1/ (D1+CCM14+M2)*=(M14M2)/4*D2)) 


H=INT(A1*286): V= INT(T3%500) 

PLOT ((H-H1),(V-V1) ,HEX(FB)) 

Hl =H: Vl =V 

NEXT Al 

EfoT (-H,-V,) 

NEXT A 

SOTO 210 

BP ww ew eee ees SUR OUR ET OPREON 34........-..+-+ 
REM ; THIS SUBROUTINE PLOTS T3 VS Al (CONST A), PARAM. L2 


INPUT "IMPEDANCE RATIO (R2/R1)=",A 

SELECT PRINT 213(64) 

PRINT " IMPEDANCE RATIO (R2/R1)=",A 

PRINT " Pole evs fictee9 DEG), PARAM L2()-0.5)" 
Beor (10) (0,-50,"+") 


mpere (10)(50,0,'"+") 

Beoe (-500,0,) 

meine = 0 TO .5 STEP .1 

meme V=O0: HlL=0: V1i=0 

mereal = 0 TO 1.56 STEP .O4 

L1 = 1/(2*SQR(1-A®A*SIN(A1) *®SIN(CA1) )) 


7S, 





960 A2 
a0 M1 
sew M2 
go B2 


1000 
ECT 0 
ILO) 748 
EOS 0 
LO4O 
iO 0 
1060 
nO) 0 
1080 
nee 
0 
0 
mez 0 
ies 0 
1140 
a 0 
eles 0 
ci? 0 
Hales 0 
ee oi 
700 
eZ AO 
ZZ 0 
es 
1240 
ikZ > 0 
acne) 
i270 
7S 0 
Zo 0 
es 
ir 0 
ina 2,0 
Levelt 
1340 
ests 0 
1360 
Leyae 
36 0 
ips 30 
1400 
1410 
1420 
1430 
L440 
1450 
1460 


ARCSIN(A*SIN(A1L) ) 
A®COS(A1)/COS(A2) 
(1/A)*COS(A2)/COS (Al) 
2*#PI*L2*L1*COS (A2) 
COS(B2)*COS(B2) 
SIN(B2)*SIN(B2) 
T3 = 1/(D1+( (M1+M2)*(M1+M2)/4#D2) ) 
H = INT(AlL"286): V = INT(T3*500) 
mor ((H=-H1), (V¥-V1),HEX(FB)) 
Hl =H: Vl=V 
NEXT Al 
Emon (=H, -V.,) 
MexT L2 
SOTO 210 
PU ve oe es ces SUBROUTINE OPTION 4 
REM ; THIS SUBROUTINE PLOTS T3 VS L (CONST A), PARAM Al 
INPUT "IMPEDANCE RATIO (R2/R10=", A 
SELECT PRINT 213(64) 
PRINT " IMPEDENCE RATIO(R2/R1)=",A 
PRINT " jia@tei) Vo L2°0-1). PARAMETER A1(22 DEG INC)" 
SELECT PRINT 413 
emer (10) (0,-50,"+") 
PLOT (10) (50,0,'"+") 
eet (-500,0,) 
FOR Al = 0 TO 1.56 STEP 
H=0: V=0: Hl=0: V1=0 
Ll = 1/ (2*#SQR(1-A*A*SIN(A1L)#SIN(A1))) 
mee 2 = 0 TO 1. STEP .02 
ARCSIN(A*SIN(A1)) 
A*COS(A1)/COS(A2) 
Gyo Gost 2 7 cos (41) 
22#PT#L2*L1*COS(A2) 
COS (B2)*COS(B2) 
SIN(B2)*SIN(B2) 
1/(D1+((M1#M2) = (M14+M2)/4"#D2)) 
H = INT(L2*500): V = INT(T3*500) 
Bior ((H-H1),(V-V1), HEX(FB)) 
Hl = H: Vl = V 
NEXT L2 
piper (-H,-V,) 
NEXT Al 
GOTO 210 
PE cee es SOCK OU E MOO TION Ss,40, ca... 2.-0ses- 
REM : THIS SUBROUTINE PLOTS T3 VS A (CONST L2), PARAM Al 
INPUT "BUBBLE LAYER/WAVELENGTH RATIO=",L 
SELECT PRINT 213(64) 
PRINT " BUBBLE LAYER/WAVELENGTH RATIO=",L 
PRINT " iieceevs ACen). PARAM Al1(22 DEG INC) 
Biov (10) (0,-50,"+") 
PLOT (10) (50,0,"+") 
PLOT (-500,0,) 


D1= 
D2 = 


eo e« eee #8 ee 8 @ @® @ @® @© © @ @ 


iso 


CO 
NO 
oak a aaa 
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7 0 
1480 
1490 
ILS he 
fe 10 
ir 2 0 
ILS eye 
1540 
BS 2 0 
iL svowe 
ifo)/ 0 
iso. S 0 
ifs 9 0 
1600 
ie 10 
1620 
iho 30 
im 00 
ey LO 
720 
i 30 
1740 
i, >0 
1760 
ee 
i SO 
i 30 
1800 
Pe 0 
38.20 
330 
1840 
1850 
ike 6 0 
ire7 0 
1880 
perc 
SoC 
Jkspae 
20 
3 0 
1940 
oo 0 
iS 6 0 
Sate 
re 0 
mgs) 0 


Oi. ae 


ee: 


FOR A= .1l 


Hex T 
al) I 
Nix £ 
GOTO 
REM 


(melowl. somo lEP .39 
Rea O0 Vy 1=0 
Oe. 3 sores 
ARCSIN(CA®SIN(CA1) ) 
A*COS(A1)/COS(A2) 
GWA) 605 (A2)/ COSCAl) 
2*#PI*L*COS(A2) 
COS CBZ) *CESt(B2) 
SINCB2)=SIN(CB2) 
1/7 (D1+¢€ (M14M2) = (M1I4M2)/4*D2)) 
mit CASO). Y = INT €@3*500) 
GCCh=hi) (Vs 1) sex CFB) ) 
ele 
A 
(-H,-V,) 
Al 
ZC 


* 
3 oe © @ © @ @ @ @ @® @® # 8 @ 


(e= 
V=0: 
= -O4 


«oe @# @ @ @® @# @® @ @® #* @® #* @&® @ @ 


meGECr R 


REM 


- THIS SUBROUTINE PLOTS T3 VS A (CONST Al), PARAM L2 


mPUT' ANGLE OF INCIDENCE=",Al1 
oe@eeCl PRINT 213(64) 
PRINT "ANGLE OF INCIDENCE=",Al1 


ERENT 'T3(0-1) VS AC0-1), 


Poke eoerbegi=.9, <1 STEPS)” 


Ramee C10) {05-50 ,"+"') 

meer (10) (50,0,"+") 

POT (-500,0,) 

eee i2 = 0 TO .5 STEP .1l 

mero =O: V=O0: VL=0 

mora = .1 TO .9 STEP .O4 

me=  1/(2*SOR(1-A*A*SINCAL)*SIN(A1) )) 
P2e= ARCSIN(CA*“SIN(A1)) 

ieee = A*=COSCA1)/COS(CA2) 

feee= ( 1/4)"*COS(A2)/COS(A1) 

ee 21 

B2 = 2*°#PI*®L*®COS(A2) 

Ml = COS(B2)*COS(B2) 

ie = SINCB2)*SINCB2) 

T3 = 1/CD1L+C(CM14M2)* (M1L4+M2)/4"D2)) 
Per rCAs500): VY = INT(CT3*500) 
mmol ((H-H1),(V=V1) ,HEX(FB)) 
mme= HH: Vl = V 

EXT A 

PLOT (-H,-V,) 

ext L2 

GOTO 210 

END 
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